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Abstract 

The  assumption  of  serial  independence  of  the  disturbance  terms  is 
the  starting  point  of  almost  all  the  work,  that  has  been  done  on 
analyzing  market  disequilibrium  models.   Under  serial  dependence  the 
usual  maximum  likelihood  estimator  (MLE)  will  be  inefficient  although 
they  may  remain  consistent.   Two  other  assumptions  about  the  distur- 
bance terras,  namely  normality  and  homoscedastici ty ,  are  also  usually 
made.   Violation  of  these  is  likely  to  make  the  MLE  to  be  inconsistent. 
In  this  paper,  we  first  derive  tests  for  serial  dependence  given  nor- 
mality and  homoscedastici ty  using  the  Lagrange  multiplier  (LM)  test 
principle.   The  likelihood  function  under  serial  dependence  is  very 
complicated  and  involves  multiple  integrals  of  dimensions  equal  to  the 
sample  size,  ruling  out  the  possibility  of  using  the  Wald  and  likeli- 
hood ratio  tests.   However,  the  test  statistic  we  obtain  through  the  LM 
principle  is  very  simple.   Next,  we  relax  the  normality  assumption 
using  the  Box-Cox  transformation  family  and  generalize  the  above  test 
to  such  non-normal  cases.   We  also  suggest  a  general  specification  test 
for  the  disequilibrium  market  model,  namely  white's  information  matrix 
test,  which  can  detect  non-normality,  heteroscedastici ty  and  serial 
dependence  jointly.   This  procedure  will  be  useful  when  there  is  not 
much  information  about  the  alternative  hypothesis.   In  the  last  part  of 
the  paper  we  present  an  analysis  of  disequilibrium  models  assuming  that 
the  disturbances  are  logistic  rather  than  normal.   This  modification 
does  simplify  the  computations  quite  a  lot  since  for  the  logistic 
distribution  a  simple  closed  form  expression  for  its  distribution  func- 
tion is  available.   The  relative  performance  of  this  distribution  com- 
pared to  the  normal  distribution  is  also  investigated. 
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1.   Introduction 

Models  of  markets  in  disequilibrium  (MMD)  have  been  studied,  both 
theoretically  and  empirically,  by  many  authors  in  recent  years  [see 
Quandt  (1982)  or  Maddala  (1983),  for  a  review].   Unlike  many  other 
econometric  models,  however,  MMD  have  yet  to  be  subjected  to  rigorous 
specification  analysis  [though  see  Quandt  (1981)  and  Lee  (1982,  1984)]. 
In  most  studies,  for  example,  disturbances  are  assumed  to  be  serially 
independent  (I),  whereas  MMD  are  designed  to  describe  time  series  data. 
Since  autocorrelation  in  economic  time  series  is  not  always  accounted 
for  by  the  exogenous  variables ,  serial  dependence  might  be  expected  to 
be  present.   Fair  and  Jaffee  (1972)  incorporated  serial  correlation  in 
their  model,  but  for  estimation  they  had  to  make  an  assumption  that 
there  was  an  equilibrium  (that  is,  both  the  demand  and  supply  were 
equal  to  the  actual  quantity  transacted)  before  every  switching  point. 
Fair  and  Kelejian  (1974)  noted  that  when  the  disturbances  are  serially 
correlated,  the  coefficients  in  the  demand  and  supply  equations  are  not 
identified.   Ameraiya  (1974)  and  Maddala  and  Nelson  (1974)  commented 
that  serial  correlation  is  a  serious  problem  and,  when  incorporated 
into  the  model,  it  renders  maximum  likelihood  estimation  com- 
putationally intractable.   The  usual  maximum  likelihood  estimators 
(MLEs)  maximize  a  likelihood  based  on  serially  independent,  identically 
distributed,  normal  disturbances.   We  continue  to  call  these  MLEs  even 
when  one  or  more  of  the  assumptions  are  in  fact  breached  by  the  data 
to  hand.   Presence  of  serial  dependence  (I)  will  render  the  MLEs 
asymptotically  inefficient,  though  they  may  still  be  consistent. 
Failure  of  the  other  assumptions,  normality  (N)  and  homoscedasticity 
(H),  may  lead  to  inconsistency  of  the  MLEs. 
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In  the  present  paper,  we  develop  some  specification  tests,  with  an 
emphasis  on  I,  for  a  particular  type  of  disequilibrium  model,  namely 
Model  1  of  Maddala  and  Nelson  (1974).   Other  MMD,  such  as  Model  2  of 
Maddala  and  Nelson,  which  assume  knowledge  of  sample  separation  will 
be  easier  to  handle.   Throughout,  the  Lagrange  multiplier  (LM)  or 
score  test  principle  is  adopted,  partly  because  of  its  good  asymptotic 
power  properties  but  mostly  because  the  likelihood  under  I  is  vastly 
more  tractable  than  that  under  I.   As  noted  by  Maddala  and  Nelson 
(1974),  the  latter  involves  multiple  integrals  of  dimension  equal  to 
the  sample  size.   [Quandt  (1981)  gives  an  expression  under  first-order 
autoregressive  (AR(1))  disturbances,  that  involves  only  univariate 
integrals  but,  as  noted  by  Lee  (1984),  and  as  we  indicate  in  the  next 
section,  this  is  not  the  likelihood.]   Thus  Wald  and  likelihood  ratio 
tests  are  not  feasible. 

In  the  next  section  we  derive  the  LM  test  for  I  assuming  NH,  our 
analysis  being  along  the  lines  of  Robinson  et  al.  (1985)  where  a  similar 
test  for  the  Tobit  model  was  proposed.   The  class  of  possible  alter- 
natives to  I  is,  in  certain  respects,  broader  than  those  in  Lee  (1984) 
where  he  considered  only  AR(1)  processes,  but  Lee's  maintained  assump- 
tion of  contemporaneously  correlated  residuals  is  more  general  than 
ours.   In  another  paper  Lee  (1982)  also  derived  a  test  for  non-normality 
(N)  with  alternatives  belonging  to  the  bivariate  Edgeworth  series  of 
distributions.   His  test  strongly  rejected  the  normality  assumption 
when  applied  to  the  housing  data  of  Fair  and  Jaffee  (1972).   The  tests 
suggested  in  Section  2  will  not  be  valid  under  N.   Therefore,  in  Sec- 
tion 3,  we  relax  the  normality  assumption,  using  the  Box-Cox  transfor- 
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mation,  and  test  for  I  given  NH.   Subsequently,  the  possibility  of 
heteroscedasticity  (H)  is  also  introduced.   Then  usiag  White  (1982) 
information  matrix  test  principle,  we  suggest  a  general  specification 


test.   This  procedure  can  be  viewed  as  a  simultaneous  test  of  INH. 
Such  a  test  is  not  likely  to  be  very  powerful  when  only  one  or  two  of 
I,  N  or  H  are  violated,  but  it  does  seem  useful  as  a  first  step  when 
there  is  little  information  concerning  possible  departures  from  the 
classical  assumptions.   A  strategy  to  be  adopted  should  any  of  the 
tests  reject  would  be  desirable,  as  would  be  a  development  of  tests  of 
N  and  H  given  I.   Unfortunately  both  objectives  are  hampered  by  the 
difficulty,  as  previously  mentioned,  of  ML  estimation  under  I.   It  may 
be  worth  investigating  whether  at  least  consistent  estimation  of  serial 
correlation  is  possible,  for  example  by  moment  methods  or  methods  simi- 
lar to  those  of  Robinson  (1982b),  which  apply  to  the  Tobit  model. 

It  is  now  well  known  that  normality  of  the  disturbances  is  not  an 
assumption  that  can  be  taken  for  granted  in  the  context  of  MMD;  nor, 
unlike  in  the  usual  regression  model,  does  it  lead  to  neat  statistical 
properties  or  to  computational  convenience.   It  is,  therefore,  worth 
considering  alternative  distributions  which  may  be  just  as  relevant  as 
the  normal  but  make  computation  easier.   Goldfeld  and  Quandt  (1981) 
considered  a  Laplace  type  distribution  which  they  called  the  Sargan 
density.   However,  for  the  simple  disequilibrium  model  this  did  not 
lead  to  much  computational  gain,  perhaps  due  to  the  "piecewise"  nature 
of  their  likelihood  function.   In  our  view,  in  this  situation,  a 
better  alternative  is  the  logistic  distribution.   The  shape  of  the 
logistic  distribution  is  very  similar  to  the  normal,  with  a  slightly 
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thicker  tail,  and,  most  important,  it  has  a  closed  from  expression  for 
its  distribution  function.   In  Section  4,  we  study  the  likelihood 
function  assuming  a  logistic  distribution,  and  generate  data  under  the 
normal,  logistic  and  Student's  t  distributions  and  make  a  comparative 
study  of  the  performances  of  the  normal  and  logistic  distributions. 

2.   Tests  for  Serial  Dependence  under  Normality 
2.1  Likelihood  Function  with  Serial  Dependence 

We  consider  the  following  simple  disequilibrium  model 

Dt  =  Xlt3l  +  Ult'  (2,1) 

St  =  X2t32  +  U2t'  (2'2) 

Qt  =  min(Dt,St)    t-1,2,...,  (2.3) 

where  D  ,  S   and  Q  denote  respectively  the  quantity  demanded, 

supplied  and  transcated  at  period  t,  x   and  x   are  column  vectors  of 

k,  and  k„  exogenous  variables  and  u,   and  u„   are  the  disturbance 
1      2    5  It      2t 

terms  such  that  E(u.  |x,  ,  s=l,2,...)  =  0,  i  =  1,2.   We  also  assume, 

it   is 

like  Quandt  (1981),  but  unlike  Lee  (1984),  that  u   and  u   are  inde- 
pendent for  al 1  t  and  s. 

Quandt  (1981)  assumed  that  the  u   are  also  AR(1)  and  set  up  the 
function 


T 

n  h  (Q  |D  ,S  ),  (2.4) 

t=l  l      t      °      ° 


where  h  (Q  Id  ,S  )  is  the  conditional  density  of  Q   given  D  ,S  ,  and 
t  xt '  o'  o  J  t  &      o'  o 

involves  only  univariate  integrals.   Quandt  (1981)  described  (2.4)  as 
the  likelihood  function,  but  in  fact  the  latter  is  of  the  form 


-5- 


T 
h(Q1,Q2,...QT)  =  h  (Q  )  H  h  (Q  |Q    ,...,Q  ).  (2.5) 

t=2 


Although  estimates  maximizing  (2.4)  may  be  consistent,  they  are  not 
efficient,  and  indeed  the  autocorrelation  coefficients  seem  likely  to 
be  poorly  identified  by  (2.4)  because  h  (Q  Id  ,S  )  *  0  as  t  *  °°. 
Perhaps  this  explains  the  seeming  contradiction  between  the  significance 
of  the  likelihood  ratio  test  for  I  and  the  insignificance  of  the  indi- 
vidual autocorrelation  coefficients,  mentioned  in  Quandt  (1981), 
though  this  might  as  well  be  due  to  the  presence  of  other  misspecif ica- 
tion(s)  such  as  N  or  a  wrong  functional  form  in  (2.1)  and/or  (2.2). 

Let  f(D,S;9)  be  the  joint  probability  density  function  (p.d.f)  of 
the  lx2T  vector  (D,S)  =  ^ ,D„ , . . . ,D  ,S  ,S2 , . . . ,S  ) ,  which  is  a  given 
function  of  D,  S  and  an  unknown  parameter  vector  8.   The  log- 
likelihood  function  is 


T      w       1-w 
L(6)  -  log  E  J  f(D,S;6)  II  (dD  )  mt(dS  )    mt 
m   R        t=l 


(2.6) 


Q 


where 


Q  =  (Q1,Q2»---»QT).  /   =  /    •••  /   . 


R    Qx     QT 


w   =  1  if  Q   =  D  , 
mt        xt     t' 


=  0  if  Q   =  S  , 
xt    t ' 


T 
and  the  sura  is  over  all  2   choices  of  (w  ,  ,w  „,..., w  „,)  ,  that  is,  over 

ml   mz      mi 

all  possible  demand/supply  sequences.   Under  our  setup  D  and  S  are 
independent,  i.e., 
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f(D,S;6)  =  f1(D;ei)f2(S;82),  (2.7) 

i    r 

where  9'  =  (99). 

We  now  introduce  I  explicitly  into  the  model.   As  a  class  of 

alternatives  to  independence  we  take,  as  in  Robinson  (1982a),  u   to  be 

generated  by  a  stationary  Gaussian  process.   The  j-th  autocorrelation 

of  u.  ,  P..(^.),  is  a  uniquely  defined,  dif f erentiable  function  of  a 
it    ij   i 

p. -dimensional  column  vector  ^ . ,  i=1.2.   Let  R.(^.)  be  the  TxT 
*i  l  11 

Toeplitz  matrix  with  (j  , j+k)-th  element  P   (^.).   If  4>  (y;M,^)  denotes 

IK.     1  X 

the  p.d.f  of  a  T-variate  normal  distribution  with  mean  u  and  variance- 
covariance  matrix  ft,  the  log-likelihood  function  can  be  written  as 

*  L(9)  =  log  Z    /  i(D;Xi   o2R(^))  ■  *  (S;X  J>  2,oh   ( #  ) ) 
m   R 


T      w       1-w 

n  (dDj  mt(dS  )   mt|   ,  (2.8) 


t-1   C       C 


:Q 


2 
where  o.    =  V(u.  )  and  X.  is  the  Txk .  matrix  with  t-th  row  x.  ,  i=l,2. 
l       it        i  l  it 

2.2  The  Test  Statistics 

We  assume,  as  in  Robinson  et  al.  (1985),  for  a  unique  value  of 
i> . ,  which  we  take  with  no  Loss  of  generality  to  be  the  vector  of 
zeros,  that  P..(0)  =  0,  i=l,2;  j=l,2,...  .   So  a  test  of  I  can  be 
expressed  as 

H  :   ■-P.  =  0,  i  =  1,2  against  H,  :   *J  *   0,  for  some  i.       (2.9) 
o     l  °        1     1 

The  alternatives  we  shall  emphasize  are  the  p.-th  order  autoregression 
[AR(p.)]  and  the  p.-th  order  moving  average  [MA(p.)]  processes.   As  in 
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the  classical  regression  case,  the  same  LM  statistic  will  result  from 
both  AR(p.)  and  MA(p.)  alternatives.   The  statistic  falls  out  quickly 
from  the  general  form  which  we  shall  derive.   Statistics  for  various 
other  types  of  alternative  can  be  readily  obtained  from  our  general 
form. 

It  is  clear  from  the  expression  (2.8)  that  the  likelihood  function 
involves  normal  integrals  of  dimension  equal  to  the  sample  size,  and  in 
general  it  cannot  be  expressed  in  terras  of  normal  integrals  of  unit 
dimension.   Therefore,  maximization  of  (2.8)  will  be  hampered  by  the 
need  to  compute  integrals  of  very  large  dimension.   Thus  the  Wald  and 
likelihood  ratio  tests  are  not  at  all  practically  feasible.   It  can  be 
shown  that  under  H  ,  when  R  (V.)  =  R_('>»0  =  I™,  where  I   is  an  identity 
matrix  of  dimension  T,  (2.8)  reduces  to  the  log-likelihood  function 
given  in  Maddala  and  Nelson  (1974,  p.  1015).   To  see  that  and  for  our 
future  use,  we  introduce  the  following  notations: 

t  t         o         t         o  t  t  t  r  i 

9A-(Br°r82'V'  V'W'  9'  =<VV- 
W  ■*i<Qt!XitBi-°?)'      1  "  1>2- 

oo 

■        Fl(Qt)  -  )    vv*itVai)dV 


and 


F2(Qt)    -   J       WXit32'°2)dSt 

^t 


For  convenience,  we  will  drop  the  symbol  Q  and  write  ^(Q*.)  ds 


f    etc.   Now  from  (2.8) 
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we  ) 


W    .  1-W    . 

mt .  .      »        mt 


log   2   /        II    {fa(D   )f2(S   )(dD   )         (dS    ) 
m     R  t=l 


■    log   I      il    (wmtFuf2t    +    U-Wmt)fuF2t> 
m  t=l 


Hog   (If       +  f     F2t)    . 

t  =  l 


(2.10) 


Since  the  log-likelihood  function  simplifies  considerably  under  H  ,  we 
turn  to  the  LM  approach  and  develop  tests  for  I. 

Let  us  consider  the  partition  9'  =  (9  .9  )  where  we  want  to  test 

A  B 

H  :   9=0.   There  are  many  asymptotically  equivalent  forms  of  the  LM 
o     B 

test  statistic  and  we  will  use  the  following  form: 


LM  =  — 77! M 


ay 


B 


B 


where  "    indicates  the  quantities  have  been  evaluated  at  the  restricted 
MLE  of  9  say  9  and  M  satisfies 

plim  [IE  C--l^i£i}|   -  M]  =  0. 

B   B  ,y 

The  above  form  is  valid  since  under  our  setup 


E    (£&1)    - 

«0  ay,; 


(2.11) 


as  shown  in  the  subsequent  derivation, 
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To   derive   LM  we    require    Che   derivatives    (under   H  ) 


-~-  log  <j)    (D;X13    .oJr^^))   =  -i_  (d-XjS   )  'R     (D-X  S   ) , 
lj 


2a; 


}  (2.12) 


■55T7  log  V8sX2s2'ff2ia<*2))  =TI  (s-x2e2)'R    (s-x232), 

2j  2o2 


where   R. .    =   (3/3^ . . )R. (0) ,    for   j    =   l,2,...,p.. 


Denoting  P...    =   (3/3y .  .  )p      (0)    and   using   (2.12)    we  have 

1 K.  J  lj        IK 


— —   [exp   L(8    )]    --i-     I    Z      p  E   /   (D  -x'    U    )(D  -x'    U    ) 

°V.  .  o  ~    2  n      l,t-u,j  t      It    1        u      lu   1 

lj  2o      t,u=l         '  J    m     R 


1 


t*u 


w 


1-w 


•    n  {f  (d  )f0(s  )(dD  )  mn(ds  )      mn)l     . 

,       1      n      I      n  n  n 

n=l 


Now 


°°  ,                                                 w                 1-w 

J  (D  -x.    3    )f    (Djf9(S    HdD    )    mC(dS .) 

Q  titlltZt                             t                                                 t 
^t 


VE2(V   * 


z   <J>    (z    ;0,a   )dz 
t    1      t'    '    1         t 


V*uBi 


+    (l-wmc)(Qt-xu3l)    •    f   (Q   )    J 


Qrx2t32 


VZt;0,°2)dzt 


Wn,tVltf2t   +   (I-V>(VXltBl)fltF2t 
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Thus 


T 
[exp  L(«  )]  >-^r  £  £  P. 


3V1J        °     2a2  C)U=i   l.t"»J 


*  '    {wmtaifltf2t  +  °-V)(VXitBl)fnF2t} 
m 


•  (w  a*f   f   +  (i-w   )(q  -x'  3  )f   F_  } 
mu  1  lu  2u       mu   u   lu  1   lu  2u 


n    {w  F.,  f_.  +  (1-w  ,  )f   F0,  1. 
i  ^t     mk  Ik  2k       mk  Ik  2k 


Following  our  derivation  of  (2.10),  we  can  write 

i     T 


"rrr exp  L(  8  )   -  -~       L     I     'p         Y   Y. 

d^T  •         o/J    .2        l,t-u,j  It  lu 
lj  2o      t,u=l       ,J 


"        tFHcf2k   +  flkF2k}> 

k=l 

k*t,u 


where    rjt    =    ^lt[0\i2t   +   (VXltVF2t]#      Therefore, 


£    £      P,  Y       Y  II      (F,,f      +f       F      ) 

.      l,t-u,j Tltrlu     .     .        Ik1 2k      lk*2k; 
t,u=l  J  k=i 


i_    [L(8    )]    --3      _^u  Jc*t^u 


3ip  o  2  T 


T-l  T 

— r-       2     P  L        V        V 

~2   .     .    lki  .     .    It    l,t-k, 

o      k=l        J  t=k+l             '         ' 
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where 


it  °lflt£2t  +   (Vxit8l)fltF2t 


1C        FU£2t   +  fltF2t 


Fltf2t  +  f!tF2t 


Under  H  ,  the  v   are  independent,  and  noting  that  the  denominator 
of  v    is  the  p.d.f  of  Q  ,  say  h(Q  ),  we  have 


«vlt)  "  / 


It 


h(Qt)    h<Vd(*t  ' 


Since  -l-Cf^F^)  =  -f^f^  „!  (Q^^B^f^,  we  can  write 


2  r   d(fltF2t}         2 

E(Vi  )  =  -oL   /  xc     dQ   =  -a  [f   F   1 

It'     1  J      dQ     ^t     1L  It  2ti« 


=  0. 


Similarly, 


2j  a^  k=l      t=k+l 


where  v 


2t 


a2fltf2t  +  (VX2t32)f2tFlt 
Fltf2t  +  fltF2t 


with  E(v2t)  =  0.   Thus 


T-l      T 

—  [L(6  )]  «-i  E  p     i      v.  v. 

3.       oJi  ik       it  i,t-k, 

i  at  k=l     t=k+l 

l 


where 


Pik 


Pikl 


ik2 


ikp  . 
l 


,    i-1,2. 
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Next 


a2      i  T"1 


-  me  >]  --  I  Plk    *   <Vir  \*±>  +  c-sg-  W-**- 

o.    k=l  t=k+l  1  l  ' 


33    3^ 
i      i 


2  T-l  T 

[L(8    )1    = -^—  L(6rt)  +  -i      Z      p             Z       {v      (—  v             ) 

^2      '    l    v    oyj  2    3v.      v    (T  2  .    ,      ik       ,       *     t\2      i.t-k' 

3a. 3^.  a.        i  a.    k=l             t=k+l             3a.         * 

li  l  l                                                   l 


i 


and  these  have  zero  expectations,  by  the  independence  of  v   over  t 

and  E(v.  )  =  0.   Of  course 
it 

E{ ^r  [L(8  )]}  =  0,  and  El    *    ;  [L(6  )]}  =  0, 

33.3^.  3az3\|;.      ° 

J   i  J   i 


for  i  *   j ,  verifying  equation  (2.11)  and 


E(— —  ,[L(8  )]}  =  0, 
1   J 


for  i  *  j.   Finally, 
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Ei-n-^-  L(d    )  -tA-  L(t>    )} 

a^.  o    0^  o 


,  T-l    T-l  T  T 

E[    ^        Epp  I  E     V     v  tv     v  ] 


a. 4        ,     .        ,      ikj    ImA   .    .  . .        ...    it   i,t-k  iu  i,u-m' 

l  k=l    m=l  J  t=k+l    u=k+l 


T-l  T 

1  9      9 

E      p.,  .p.,  .    E[      ^      vf   vf  ] 


a. 4  Ikj    iki      \  .  _,    It   i,t-k' 

l      k=l  J  t=k+l 


T-l  T 

Z      P...P...        2      A      A  where    A        =  -t  E(v2   ).  (2.13) 

k=1      ikj    ik*   t=k+1    it    i,t-k  it        Q2  it/ 

i 


2 
A  closed  form  expression  for  E(v   )  does  not  appear  to  be 

available,  however,  so  we  replace  (2.13)  by 

T-l  T 

V  =  /,pikjpiwcik>  cik =a-    I ,  vitui,t-k> 

J    k=l   J  a. 4  t=k+l 


2  "2 

where  v.   is  v.   with  tf.  and  0.  replaced  by  their  MLEs  3.  and  a. 

it     it       l      l    K  J             11 

The  LM  statistic  is  thus  given  by 


2   , 

LM  =  l    l.H.l. , 
i-1  2  J  * 


,   T-l  T    -   - 

where  L.    =- —  £   P.,d.,  with  d.,  =   Z        *>  v      and  M.  =  ((M.  .,,)). 

1    a. 2  k=l   lk  lk  lk   t-k+1   1C  1,t_k      *      1J" 

l 

When  the  null  hypothesis  H   is  true,  under  certain  regularity 
condi  tions 

V 

LM  +  x2       •  (2.14) 

Pl  +  P2 
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These  regularity  conditions  may  be  established  using  the  results  of 
Hartley  and  Mallela  (1977)  and  combining  the  conditions  of  Robinson 
et  al.  (1985)  on  the  P..(>K). 

We  now  deduce  the  form  of  LM  for  testing  against  AR(p.)  and  MA(p.) 
alternatives.   For  both  the  AR(p  )  and  MA(p.) 

P.,  .  =1  for  k  =  j 
lkj  J 

=  0  for  k  *  j, 


so  that  in  either  case,  the  LM  test  statistic  reduces  to 

LM(p  p  )  =  T  I        I      rf  (2.15) 

1      l  i=l  k=l   llc 

1   dik 
where  r.,  =  —  -s- .   Note  that  (2.15)  is  analogous  to  the  familiar 

lk    ,-   2  l/n 

/T  c  cf/2 

i  ik 

Box-Pierce  portmanteau  statistic  for  testing  for  serial  correlation. 


2.3  Empirical  Illustration 

We  apply  the  test  to  the  housing  start  model  of  Fair  and  Jaffee 
(1972).   Their  versions  of  (2.1 )— (2.3)  are  [see  Maddala  and  Nelson 
(1974,  p.  1025)] 


Dt    =a0+Vlt   +U2Z2t    +V3t   +   Ult'  (2'16) 

St   =   30   +  Vlt   +  d2Z4t   +   y3Z5t   +   Vet    +   U2f  (2'17) 

Qt    =   min(D    ,S    ),  (2.18) 


-15- 


where  Q   is  the  observed  number  of  housing  starts  in  period  t,  z   is 

a  time  trend,  z   is  a  measure  of  the  stock  of  houses,  z~   is  the 

mortgage  rate  lagged  two  periods,  z,   is  the  six  month  moving  average 

of  the  flow  of  private  deposits  in  savings  and  loan  associations 

(SLAs)  and  mutual  savings  banks,  lagged  one  period,  z   is  the  flow  of 

borrowings  by  the  SLAs  from  the  Federal  Home  Loan  Bank  lagged  two 

periods,  and  z   =  z  .   Using  the  Fair  and  Jaffee  (1972)  data  with 

o t    J , t+1 

a  sample  of  size  126,  Maddala  and  Nelson  (1974,  Table  II)  obtained  the 
MLEs  which  are  reproduced  below: 

n      -    I    an     ai      ao     aT     Sn     t*i     3o    S-j    *,             a?    ao 
Parameter     0 1 2 3 0 12    3    4 1    2 

Estimate  j  223.74  2.52  -.022  -.090  15.55  -.153  .053  .053  .093  350.0  80.2 


Lee  (1982)  found  evidence  of  no  contemporaneous  correlation  between  the 

u.   for  these  data,  so  the  setting  for  an  LM  test  for  I  seems  appropriate 
here.   We  use  the  above  estimates  and  obtain  the  following  components 
of  the  LM  statistic  given  in  (2.15): 

Trll  =  4*4733       Tr21  =  23.9101      It1  =  6.0109 

Tr^   =   .4636       TrJ_  =  10.3031      TrJ,  =   .9269 
i~l                                            11                                           26 

Tr^  =  1.7940       Tr^  =   9.0344      Tr^  -   .3906 

Tr^4  =   .0003       Tr^  =   6.7310      Tr^g  =   .2724 

It  seems,  the  length  of  the  disturbance  AR  process  is  short  for  the 
demand  equation  whereas  for  the  supply  equation  it  is  quite  the  oppo- 
site.  Using  the  above  values  we  can  form  a  number  of  test  statistics, 
e.g.  , 
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LM(1,0)  =  4.4733  LM(0,1)  =  23.9101 
LM(1,1)  =  28.3834  LM(1,2)  =  38.6865 
LM(2,1)  =  28.8470        LM(2,2)  =  39.1501 

2 
The  X  critical  values  for  1,  2,  3,  and  4  degrees  of  freedom  at  one  per- 
cent significant  level  are  respectively  6.63,  9.21,  11.34,  and  13.28. 
Therefore,  we  strongly  reject  the  null  hypothesis  of  no  serial 
dependence.   The  inference  remains  the  same  even  if  we  consider  higher 

values  of  p  and  p?.   Here  we  should  note  that  our  sample  size  is 

2 
126,  and  the  asymptotic  X   critical  values  are  unlikely  to  be  appro- 
priate.  However,  for  this  particular  example,  the  statistic  values  for 
the  supply  equation  are  so  large  that  the  conclusion  of  u   being 
serially  correlated  seems  justified  insofar  as  the  maintained  assumption 
of  correct  specification,  normality  and  horaoscedasticity  are  valid. 

3.   Tests  for  Serial  Dependence  in  the  Presence  of  Non-normality  and 
A  General  Specification  Test 

3.1  The  Test  Statistics  for  Serial  Dependence  under  Non-normality 

In  the  last  section,  a  test  was  derived  under  the  assumption  that 

the  disturbances  are  normally  distributed.   When  this  assumption  is 

violated,  the  usual  MLEs  will  be  inconsistent  and  the  test  will 

be  invalid.   We  nowallow  for  N  using  the  Box-Cox  transformation 

family,  generalizing  our  previous  test. 

We  consider  the  model  as  given  in  (2.1)— (2.3)  except  now  assume 

that  initially  the  u.  's  are  not  normally  distributed  but  a  Box-Cox 

it  ' 

transformation  of  the  dependent  variables  restores  normality. 
Therefore,  now  we  have 


where 
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Dt    "  Vl  +  Ult'  (3'1} 

(V 

St    =X2t62+U2t'  (3-2) 

Qt  =  mln(Dt,St),  (3.3) 


(*  )  D /  -  1 

D,    =  T(Dt,A  )  --S^-   if  \   *  0, 


(3.4) 


=  log  Dt    if  X±   =  0, 

(V        (V         2 
and  the  p.d.f.  of  D     is  <J>  (D    ;  X^.,  a  (R  (♦- )).   Similarly  we 
(A  )  i         i  i   1   i   1 

define  S     and  its  p.d.f. 

Strictly  speaking,  the  u   ' s  cannot  be  normal  because  of  the 

nature  of  the  transformation  in  (3.1)  and  (3.2).   For  (3.4)  to  be  well 

defined,  D  and  S   should  always  be  positive.   This  implies  the  u.  's 
'   t      t  it 

cannot  be  normal.   An  appropriate  assumption  in  this  situation  would 
be  that  the  u   's  are  truncated  normal  [see  Poirier  (1978)] 

T(0,A.)  -  x!  8.  <  u.   <  T(°°,A.)  -  x!  3.  . 
l      it  l     it       '  l      it  l 

However,  here  we  are  dealing  with  macro  demand  and  supply  functions; 

t 

therefore,  the  x.  tf.'s  are  positive  and  can  be  expected  to  be  much 
'       it  l 

larger  than  ju.  |.   Hence,  for  all  practical  purposes,  we  can  neglect 


the  effects  of  truncation. 

=  (  l5  .  Q 

A 


d.    -  Wv    ar  A1>  t$2,  a2,  X2),  ag  =  (V;L,  *2), 


d'  =  (^a»  ^n)  and  the  log-likelihood  function  for  the  model  (3.1)— (3-3) 

A   B  ° 

can  be  written  as 
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L(8)    =   log   2    /   CT(D;    X8,    Ap    a*R  (ij;   )) 
m      R 


CT(S;    X232,    X2,    o^R2(^2)) 


w  1-w  4 

mty  ._    v        mt 


n  (dDt)  mu(dst) 


where 


€T(D;    Vl'    V   aiW}    =   (27r)'T/2|a2Ri(^)| 


i      <V  -i   (V  T  V1 

•    exp    [--±j  (D  "Vl),Il(*l)      (D  ~Vl)]    '  Dt        " 


2a' 


t=l 


and  similarly  we  define  ^T(S;  X232'  X    a  R  ( ip  ) ) . 
Now  let 


9'.  =  (^!,  a2,  X.,  ij/Y  d'.   =  (g'   a2,  X.,  0»),  i  =  1,2, 
1     1    1 '   1    1     10     1    1    1 


1    » 


e-  -  (e   82),  and  eo  -  e   ,  > 


Then  proceeding  as  in  the  previous  section,  we  have 

^ «- «v'i -i  ti -i.t,.j ;  »tv  ■ 

1 


x  ,rf,)(D      J      -    x      3    ) 
1 1    1         u  lu    1 


t*u 


w 


mn, 


1-w 


II       (f.(D    )f9(S    )(dD    )         (dS    ) 
1      n      L      n  n  n 

n=l 


mn 


(V     ■       2   V1 

where    fjCD^    =   ^^         ;    x}  ^  ya  1  )Dt         .      Now 
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V^1' -  xit8i)f1(vf2(st)(dDt)Wmt(dst)1  v| 

Qt  lQc 


where    f        -  ♦    (Q      X    ;    «J   »      1)    and   I       -  /  ♦  (.   j    «jt»2.«*>d.   . 

since   /      (D      2    -   *ltVfl(Dt)dDt   =  Vl(Qt         ;    xltSl,1)# 


Therefore,  the  derivations  and  the  resulting  test  statistic  are 

essentially  the  same  as  in  the  previous  section  with  new  definitions 

of  f,(Q«.)t  F-(CL)  and  v,  .  i=l,2.   For  this  case  v   will  be 
i  xt    i  xt       it  It 

2-          (V 
<Cf.  f.  +  (Q  -1   -  x'  3,)f,  F0 
1  It  2t    vt It  1   It  2t 

It  F   f    +  f   F  ' 

It  2t    12  2t 


with  a  corresponding  expression  from  v   .   We  can  also  show  that 

E(v.  )  =  0,  and  clearly  the  v.   are  independent  under  I. 
it  3  it         v 

Under  the  above  framework  we  can  easily  introduce  H,  and  test  for 
I  and  H  given  N  or  test  for  I,  N  and  H  simultaneously.   One  way  to 

<V 

incorporate  H  would  be  to  assume  that  the  p.d.f.  of  D     is 
(A  ) 

^T(D  l    ;X1B1,«1)  where  Uj  =  ^|/2R1(^1)^]/2  where  ^  =  diag(^  +  z^) 


with  z   and  ol      being  known  and  unknown  vectors.   Similarly,  we  can 
incorporate  H  in  the  distribution  of  S     .   Under  this  general  set- 
up, for  a  simultaneous  test  the  null  hypothesis  would  be  H  : 

*.  =1,  a   =0  and  i|» .  =0,  i=1.2.   Following  the  above  lines,  the  LM 

11  l  ° 

test  statistic  can  be  easily  derived  and  this  can  be  viewed  as  a 
general  model  specification  test.   However,  if  we  are  interested  in  a 
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simpler  overall  model  specification  test  without  any  specific  class  of 
alternatives  in  mind,  it  seems  a  wiser  approach  would  be  to  use  the  White 
(1982)  information  matrix  test.   In  the  next  subsection  we  outline  the 
test  procedure. 


3.2  A  General  Model  Specification  Text 

Let  us  start  with  the  model  (2.1)-(2.3)  under  the  assumption  that 

2 
the  u   are  iid  N(0,o.),  i*l,2.   We  are  interested  in  testing  the 

'   2   '   2 

validity  of  this  model.   Let  us  denote  01  =  (3,,  a   3   a  )  and  let 

the  dimension  of  9  be  p.   Under  this  set-up  the  density  function  for 
the  t-th  observation  is  given  by  [see  equation  (2.10)] 


g(n  • 6)  =  p  f   +  f  F 
8VV  ;    It  2t    ltr2t  * 

If  the  model  is  correctly  specified  then 

32log  g(Q;B)1        j  log  g(Q;0)   ,3  log  g(Q;9) 
^g  [~  3  0  30'    J     g  K     30       ;  v 30       ;  J' 


where  E  denotes  the  expectation  assuming  g( • , " )  is  the  true  probability 
density  function.   White's  information  matrix  test  essentially  tests 
the  above  identity.   Such  a  test  can  be  based  on  the  following  quan- 
tities [see  White  (1982,  p.  9)] 

r±  *  *  * 

3  log  g(Qt;«)   <Uog  g(Qt;0)   dlog  g(Qt;6) 

wtz  =  5TT5 + IB 9  a3 »        (3'6) 

j   k  j  k 

j    k  =  1,  2,  ...,  p,   £  »  1,  2,  ...,  p(p+l)/2, 
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where  9  and  9.  are  respectively  the  MLE  and  j-th  element  of  6.   Asymp- 
totic moments  of  the  quantities  in  (3.4)  under  the  null  hypothesis  are 
quite  complicated.   However,  Chesher  (1983)  provides  a  simpler  way  to 
calculate  an  asymptotically  equivalent  version  of  the  test  statistic. 
Let  W  be  a  Txq  matrix  as  defined  in  (3.4)  with  q  =  p(p+l)/2.   Then  the 
test  statistic  can  be  calculated  as 

S  =  1_'  W(W»W)~1W,1_  -  T  *  R?  .  w, 

2 
where  J_  is  a  Txl  vector  of  ones  and  R,  .    is  the  uncentered  coef- 
ficient of  detemination  from  the  regression  of  a  vector  of  ones  on  W. 

2 
Under  the  null  hypothesis,  S  is  asymptotically  distributed  as  X  with 

q  degrees  of  freedom.   To  calculate  S  all  we  need  are  the  first 

and  second  derivatives  of  g(Q  ;8)  evaluated  at  8  =  9 ,  and  these  are 

given  in  Maddala  and  Nelson  (1974,  pp.   1016-1018).   Therefore,  it  is 

not  difficult  to  carry  out  this  test.   Here  we  should  mention  two  points. 

2 
First,  the  X   degrees  of  freedom  will  be  quite  large  even  for  a  small 

model  so  the  test  may  not  be  very  powerful.   One  possible  solution  is 

2 
to  consider  only  elements  with  j=k  in  (3.6).   Then  we  will  have  a  x 

test  with  p  degrees  of  freedom.   Second,  this  is  a  pure  significance 

test  and  therefore,  once  the  null  hypothesis  of  no  misspecif ication  is 

rejected,  it  does  not  provide  any  information  regarding  the  direction 

of  departure  from  the  tested  model.   One  possible  strategy  would  be  to 

apply  one-directional  tests  such  as  the  Lee  (1982)  non-normality  test 

and  a  test  for  serial  dependence  separately  with  low  significance 

levels. 
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4.   Disequilibrium  Analysis  with  Logistic  Distribution 
4.1  Preliminaries 

One  problem  with  the  log-likelihood  function  based  on  the  normal 
distribution  is  that  its  evaluation  requires  computation  of  normal 
integrals  for  which  no  closed  form  expressions  are  available.   Also, 
earlier  experience  in  estimating  disequilibrium  models  has  revealed 
difficulties  in  obtaining  convergence  using  a  Gauss-Newton  type 
algorithm.   Therefore,  there  is  a  need  to  look  for  alternative  distri- 
butions which  will  lead  to  well-behaved  and  computationally  simpler 
log-likelihood  functions.   In  this  section,  we  will  explore  the  possi- 
bility of  using  the  logistic  distribution  in  analyzing  disequilibrium 
models. 

If  a  random  variable  z  follows  a  univariate  logistic  distribution, 
its  p.d.f  h  and  distribution  function  H  can  be  written  as 

,     exp  {-(r^-)} 
h(z;u,a)  =  j- 


[1    +  exp  {-(^))]2 


and 

H(z;m,o)  = 


1  +  exp  i-(V-)} 


-00<z<0°,  -co<m<qo,  0  <  a  <  <", 

2  *2 
respectively.   It  is  well  known  that  E(z)  =  M,  V(z)  =  cf  — —  t 
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the  coefficient  of  skewness  -  0.0  and  Che  coefficient  of  excess  =  1.2. 
Therefore,  the  logistic  distribution  is  very  much  like  the  normal  but 
with  a  slightly  thicker  tail. 

4.2  Log-likelihood  Function  and  Estimation 

If  we  assume  that  D  and  S   follow  two  independent  logistic 
distributions  then  the  log-likelihood  function  for  the  disequilibrium 
model  can  be  written  as 


L(6)    =      I   log   G 
t=l 


where 


G      = 


exp    l-(Qt     '  Xlt3l)}  exp{-(Qt   "  X2ttf2)} 

1  1  2 


C        °1 Qt     -Xit3lul2        n    +  .    ^t-X2t32 

)}]  [1   +  exp    {-( — 

1  2 


r    ,Q      -   x'    3 .  x ,  ,    .0      -   x'    doN . 

exp    l-(xt           It    1)}  exp{-(xt          2t   2)) 

a  a 

1  1                                        2 


Q      -   x*    8  a„  Q      -   x'    ci 

n  r    /    t  It    lu,  2      .  ,    ,t  2t    2,  ,  ,2 

[1   +   exp    {-( )}]  [1    +   exp    {-( )}] 

1  2 


Q  -  x!  3  . 
t    i  t  i 

By  denoting  h.   =  a.   =  exp  (-h.  ),  i  =  1,  2 ,  we  have 

it       a.        it     *    ,      It 


L(6)  =  I    [log  I    (a.(i+a.  )}  -   I  {h.   +  21og(l+a. J  +  log  o.}}         (4.1) 
t=l      1-1   x     1C      1-1   1C  ll 
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The  above  log-likelihood  function  is  much  simpler  than  (2.10),  or  the 
expression  given  in  Goldfeld  and  Quandt  (1981,  p.  149)  based  on  the 
Sargan  density.   It  is  clear  that  L(9)  given  in  (4.1)  is  continuously 
dif f erentiable.   By  choosing  appropriate  parameter  values  such  as 
a   =  0,  x  3   =  Q  ,  o     *   0,  x  3  *   Q  for  t  *   1,  it  is  seen  that 
L(9)  is  unbounded.   This  property  also  holds  for  the  log-likelihood 
function  based  on  the  normal  distribution  [see  Quandt  (1982,  p.  13)]. 
It  would  be  interesting  to  investigate  the  performance  of  the  two 
log-likelihoods  given  in  (2.10)  and  (4.1)  in  terms  of  computation  and 
statistical  properties  of  the  estimates  when  the  data  are  generated 
from  different  distributions.   This  is  what  we  set  out  to  do  in  the 
next  section  through  simulation  experiment. 

4.3  Simulation  Results 

The  simulation  experiment  was  carried  out  using  the  housing  start 
model  given  in  (2. 16)-(2. 18)  and  the  Fair  and  Jaffee  data  set  on  Z 
variables.   The  disturbance  terras  u  and  u  were  generated  as  the  nor- 
mal, logistic  and  Student's  t  with  7  degrees  of  freedom.   For  the 

population  parameter  values  we  used  the  MLEs  of  Maddala  and  Nelson 

2  2 

(1974)  except  we  set  ot   =  171.8,  o      =   40  and  o      =  30.   These  changes 

were  made  to  make  the  task  of  maximizing  the  log-likelihood  function  a 
little  easier.   Observations  on  D  and  S  were  obtained  using  (2.16) 
and  (2.17),  and  the  (2.18)  provided  the  Q   values.   We  generated  100 
samples  each  of  size  126  under  the  above  three  distributions  and  esti- 
mated the  model  assuming  that  the  underlying  distribution  is  either 
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normal  or  logistic.   We  employed  a  Gauss-Newton  type  iteration  scheme 
to  obtain  the  MLEs ,  namely, 


(r)  =  Vl)  +  6<A'A>"lA,I 


,  r  >  2, 
(r-1) 


where  9   v  is  the  value  of  9  at  the  r-th  iteration,  S  is  a  scalar 
(r) 

chosen  at  each  iteration  step  and  A  =  ((A  ))  =  (( oL  ( 9)/ 39 . ) ) ,  L  (9) 
being  the  log-density  function  for  the  t-th  observation.   We  experi- 
mented with  many  starting  values  9    for  the  iterations,  such  as 
estimates  obtained  from  the  Directional  Method  I  considered  by  Fair  and 
Jaffee  (1972)  and  final  iterated  values  from  an  arbitrary  replication. 
(In  many  cases,  the  iterations  converged  to  different  points  which  were 
close  to  the  starting  values  indicating  the  presence  of  locally 
optimal  points).   After  much  trial  and  error,  we  found  that  the  latter 
option  gave  better  results.   For  the  logistic  log-likelihood  we  used 
numerical  derivatives,  whereas  for  the  normal  log-likelihood  we  had  to 
provide  the  analytical  derivatives  in  order  for  the  iteration  process 
to  converge.   This  indicates  one  advantage  of  using  the  logistic 
distribution. 

After  obtaining  the  estimates  from  100  replications  we  calculated 
the  root  mean  square  errors  (RMSE's)  for  different  cases,  and  these 
are  reported  in  Table  4.1.   The  results  are  quite  mixed  and  Che  RMSE's 
do  not  differ  much  for  the  six  cases.   Use  of  the  correct  distribution 
does  not  always  lead  to  lower  RMSEs .   This  may  be,  apart  from  con- 
vergence to  suboptimal  points,  due  to  finite-sample  errors  and/or 
the  small  number  of  replications.   Under  the  correct  setup,  the 
logistic  distribution  does  slightly  better  than  the  normal,  which  is 
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Table  4.1 

Comparison  of  Root  Mean  Square  Errors 

True  Distribution       Normal Logistic        Student's  t 

Assumed  Distribution    Normal  Logistic  Normal  Logistic  Normal  Logistic 

Parameters  a         10.6150   9.5568  10.6073  10.2584  10.5496   11.3700 


al 

.1532 

.1382 

.1531 

.1487 

.1523 

.1654 

a2 

.0014 

.0013 

.0014 

.0014 

.0014 

.0015 

a3 

.1750 

.1768 

.1750 

.1770 

.1751 

.1774 

6o 

.9774 

.8912 

.9768 

.9613 

.9715 

1.0711 

h 

.7816 

.6293 

.7804 

.6767 

.7756 

.6634 

h 

.2860 

.2224 

.2856 

.2217 

.2839 

.2186 

33 

.2558 

.1950 

.2553 

.1951 

.2536 

.1941 

\ 

.3521 

.3419 

.3507 

.3503 

.3481 

.4123 
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not  the  case  when  we  use  the  normal  distribution  correctly.   When  the 
distribution  is  misspecif ied,  the  performance  of  the  normal  distribution 
is  better,  though  the  differences  are  not  substantial.   Therefore, 
bearing  in  mind  its  computational  advantages,  the  logistic  distribution 
might  be  preferred  to  the  normal  in  the  analysis  of  market  disequilibrium 
models  with  independent  disturbances. 


-28- 


Ref erences 

Chesher,  A.  (1983),  "The  Information  Matrix  Test,  Simplified  Calcula- 
tion Via  A  Score  Test  Interpretation,"  Economics  Letters,  13, 
45-48. 

Fair,  R.  C.  and  D.  M.  Jaffee  (1972),  "Methods  of  Estimation  for 
Markets  in  Disequilibrium,"  Econometrica,  40,  497-514. 

Fair,  R.  C.  and  H.  H.  Kelejian  (1974),  "Methods  of  Estimation  for 
Markets  in  Disequilibrium:   A  Further  Study,"  Econometrica,  42, 
177-190. 

Goldfeld,  S.  M.  and  R.  E.  Quandt  (1981),  "Econometric  Modelling  with 
Non-normal  Disturbances,"  Journal  of  Econometrics,  17,  141-155. 

Hartley,  M.  J.  and  P.  Mallela  (1977),  "The  Asymptotic  Properties  of 
the  Maximum  Likelihood  Estimator  for  a  Model  in  Disequilibrium," 
Econometrica,  45,  1205-1220. 

Lee,  L-F.  (1982),  "Test  for  Normality  in  the  Econometric  Dis- 
equilibrium Market  Model,"  Journal  of  Econometrics,  19,  109-123. 

Lee,  L-F.  (1984),  "The  Likelihood  Function  and  A  Test  for  Serial 

Correlation  in  A  Disequilibrium  Market  Model,"  Economics  Letter, 
14,  195-200. 

Maddala,  G.  S.  (1983),  Limited-Dependent  and  Qualitative  Variables  in 
Econometrics ,  Cambridge:   Cambridge  University  Press. 

Maddala,  G.  S.  and  F.  D.  Nelson  (1974),  "Maximum  Likelihood  Methods 
for  Models  of  Markets  in  Disequilibrium,"  Econometrica,  42, 
1013-1030. 

Poirier,  D.  J.  (1978),  "The  Use  of  Box-Cox  Transformation  in  Limited 
Dependent  Variable  Models,"  Journal  of  the  American  Statistical 
Association,  73,  284-287. 

Quandt,  R.  E.  (1981),  "Autocorrelated  Errors  in  Simple  Disequilibrium 
Models,"  Economics  Letters,  7,  55-62. 

Quandt,  R.  E.  (1982),  "Econometric  Disequilibrium  Models,"  Econometric 
Reviews ,  1,  1-63. 

Robinson,  P.  M.  (1982a),  "On  the  Asymptotic  Properties  of  Estimators  of 
Models  Containing  Limited  Dependent  Variable  Variables," 
Econometrica,  50,  27-41. 

Robinson,  P.  M.  (1982b),  "Analysis  of  Time  Series  from  Mixed  Distribu- 
tions," Annals  of  Statistics,  10,915-925. 


-29- 


Robinson,  P.  M. ,  A.  K.  Bera  and  C.  M.  Jarque  (1985),  "Tests  for  Serial 
Dependence  in  Limited  Dependent  Variable  Models,"  International 
Economic  Review  (forthcoming). 

White,  H.  (1982),  "Maximum  Likelihood  Estimation  of  Misspecif ied 
Models,"  Econometrica,  50,  1-25. 


D/311 


